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UNIVERSALITY, SCALING AND TRIVIALITY IN A 
HIERARCHICAL SCALAR FIELD THEORY 
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Using polynomial truncations of the Fourier transform of the RG transformation of Dyson's hierarchical model, 
we show that it is possible to calculate very accurately the renormalized quantities in the symmetric phase. 
Numerical results regarding the corrections to the scaling laws, triviality, hyperscaling, universality and high- 
accuracy determinations of the critical exponents are discussed. 
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1. Introduction 

Probing shorter distance physics by direct pro- 
duction experiments is becoming increasingly dif- 
ficult. In the 21-th century, we may have to look 
for small effects at accessible energies in order to 
learn about physics at higher energy. If this sce- 
nario is correct, it is important to develop meth- 
ods of calculations which can outperform the ac- 
curacy of the Monte Carlo methods. 

We present here some results showing that very 
simple algebraic methods can be used to calcu- 
late very accurately the renormalized quantities 
of Dyson's hierarchical model |Q . The RG trans- 
formation of this model is closely related 0] to 
Wilson's approximate recursion formula [pIT In 
both cases (which are examples of "hierarchical 
approximations"), the RG transformation maps 
a local potential into another local potential and 
there is no wave function renormalization or gen- 
eration of higher derivative terms. Such a simpli- 
fied version of the RG is only a qualitative ap- 
proximation of the one which holds for nearest 
neighbor interaction lattice models. However, if 
we consider the values of the critical exponents, it 
is a much better approximation than the gaussian 
model. 

As far as the zero momentum Green's functions 
are concerned, the numerical treatment that we 



proposed in Ref. j|] is tantamount to a closed 
form solution. The method advocated replaces 
the evaluation of multiple integrals (which in 
many lattice models can only be performed with 
the MC methods) by the evaluation of a single 
integral followed by purely algebraic manipula- 
tions. The computer time involved in this proce- 
dure grows only like the logarithm of the number 
of sites. These great numerical advantages lead 
us to reconsider the question of the improvment 
of the hierarchical approximation. This is a hard 
bookkeeping and computational problem which 
we plan to attack in the future. 

Dyson's model couples the main spin in boxes 
of size 2 l with a strength (|) , where c is a free pa- 
rameter set to 2( D ~ 2 '/ £> in order to approximate 
a nearest neighbor scalar model in D-dimensions. 
A more systematic presentation of the model can 
be found for instance in Ref. Models with 

fermi fields can be constructed similarly by re- 
placing D — 2 by D — 1 in c. In addition one needs 
to specify a local measure Wo (</>), for instance of 
the Landau-Ginzburg type {W {4>) = e -^ 2 - B * 4 ) 
or of the Ising type. Under a block spin transfor- 
mation, the local measure changes according to: 



W n+ i{4>) oc e*U 
x / dcj>'W n ( 
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This recursion formula can be reexpressed in 
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Fourier representation. 
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where R n {k) is the Fourier transform of W n (<p) 
with an appropriate rescaling (^]. 

2. The Calculational Method 

Recently, it was found Q that finite dimen- 
sional approximations of degree l m ax- Rn(k) = 



1 + a n .\k + a n ok 



+ a n-l ma!c k 



2/„ 



provide 



very accurate results in the symmetric phase and 
moderately accurate results in the broken phase. 
In order to fix the ideas, one can calculate the 
critical exponent 7 with three decimal points us- 
ing l max = 15 and with 13 decimal points with 

'max 

= 50. The recursion formula for the a n m 



reads a n+ \ t 
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The initial condition for the Ising measure is 
R = cos(k). For the Landau- Ginsburg measure, 
the coefficients in the k— expansion need to be 
evaluated numerically. This is the only integral 
which needs to be calculated, after we only have 
algebraic manipulations. The effects of finite di- 
mensional truncations decay faster than exponen- 
tially. If x denotes the susceptibility calculated 
with l max = I, then 0] 
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The volume effects can be controlled arbitrarily 
well in the symmetric phase where all the calcu- 
lations which follow have been performed. The 
number n(/3, A) of iterations necessary to calcu- 
late the susceptibility at fixed /3, with a relative 
precision | x " + ^~ x " | = A can easily be estimated 
as 

n(/3,A) = -(^)(Log 2 (A) +1 Log 2 (f3 c -(3)) .(5) 

The main source of error comes from the round- 
off errors. Near criticality, R n (k) spends many 



iterations near the fixed point and the errors are 
amplified in the unstable direction. A simple cal- 
culation [Q shows that if 6 is a typical round-off 
error (e.g. 1CP 16 in double precision), then the 
relative error on the susceptibility is of the order 
(f3 c — (3)~ 1 8. Numerical experiments confirm this 
estimate, however the detailed statistical prop- 
erties has some intriguing features (non-gaussian 
distributions) which remain to be understood. 

3. Numerical Results 

The method described above can be combined 
with conventional expansions. In particular it al- 
lows us to calculate several hundred coefficients 
of the high-temperature expansion. In D = 3, it 
was found |(| that if we use a parametrization of 
the form 



X = (ft* - /3r 7 (A, + A X {I3 C - /3) A + ....) 



(6) 



the A , Ai are in general log-periodic functions of 
the form J2i a i(Pc ~ fi) 7 "^ where A is the largest 
eigenvalue and ai ^ . In D = 4, the loga- 
rithmic corrections to the mean-field scaling laws 
can also be obtained from the high-temperature 
expansion. Minimizing the errors on t(m) = 
{{T m p c — l)m) _1 , where r m is the ratio of two 
successive HT coefficients, for 300 < m < 400, we 
found H that 



x ~tf e -p)-^(M\ln(Pc-0)\) p + Ai) 



(7) 



with 7 = 0.9997 and p = 0.3351. This result is in 
good agreement with the classic field-theoretical 
result 7=1 and p = 1/3. 

Eq. (|3|) can be used to calculate numerically the 
renormalized coupling constants. Using the nota- 
tion M n for the total field ^ <fi x inside blocks of 
side 2™ and < ... >„ for the average calculated 
inside these block, we define the dimensionless 
renormalized coupling 



< tf K 4 > n +3(< Ml > n f 
2 «(<^)f+2 



(8) 



The numerator scales like 2™ while its two terms 
scale like 2 2 ™ and as n increases, more and more 
significant digits get lost in the subtraction pro- 
cedure. It is nevertheless possible to stabilize B 
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several digits of A4. We found that for D = 3, 
A4 reaches a finite non-zero limit at criticality 
A| = 1.92786 for a Ising measure. This property 
is sometimes refered to as hyperscaling. When 
one approaches criticality (which can easily be re- 
formulated as "when the cut-off becomes large"), 
this limiting value is approached according to the 
approximate law 



A4 — Xt 



1.68 x [fi c - py 



-0.43 



(9) 



In D = 4, one can check the "triviality" of the 
theory For an Ising measure, A4 reaches zero 
when p tends to p c according to the approximate 

law a 



(I/A4) ~ -1.96 - 0.746 x ln(p c - p) 



(10) 



One can also use Eq.(|3|) to calculate approxi- 
mate fixed points of the RG transformation. For 
this purpose, we start with an arbitrary measure 
and we fine-tune P until R n {k) stabilizes for a 
large number of iterations. This can be monitored 
in terms of ratio of successive coefficients. The ap- 
proximate fixed point so obtained are fixed points 
for a particular value of P, however, it is possible 
to reabsorb this dependence by a rescaling of k. 
We have followed this procedure for a large class 
of models @, namely W {<j>) oc exp -(irn 2 4> 2 + g 4> 2 n 
with rri 2 — ±1 (single or double-well potentials), 
p = 2, 3 or 4 (coupling constants of positive, zero 
and negative dimensions when the cut-off is re- 
stored) and g = 10 or 0.1 (moderately large and 
small couplings). All the approximate rescaled 
fixed points R*(yfpZk) we found turned out to be 
very close to a a universal function 



U(k) = 1. - 0.35871 * k 2 + 0.05354 * k 4 



(11) 



This universal function coincide in very good ap- 
proximation with the function which can be ob- 
tained from the fixed point constructed with great 
accuracy by Koch and Wittwer S. Namely, we 



found that \Sui\ < 



5xir 

V.2> 



where ui are the co- 



efficients of U. The use of this fixed point al- 
lows a very accurate determination |10| of the 
critical exponents appearing in Eq. (|): 7 = 
1.2991407301586 andA = 0.4259468589881. Di- 
rect fits of the susceptibility confirm 11 decimal 
points of 7 and 5 of A. 



Work in progress involves the improvement of 
the hierarchical approximation, development of 
accurate methods in the broken phase and tests 
of perturbation theory. We are also considering 
hierarchical Fermi-Bose systems in order to test if 
it is possible to construct fully non-perturbative 
model without fine-tuning. In other words can 
the addition of hierarchical fcrmions with a suit- 
ably chosen set of bare couplings to bosons drive 
the boson measure "naturally" toward the stable 
manifold? Another question being investigated 
is: given the fixed point of Dyson model, can we 
calculate the susceptibility away from criticality? 
Or in other words, does the result of Koch and 
Wittwer || solve the model ? The answer to this 
question seems to be yes in low l ma x cases. 
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